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Abstract 

To solve the 'TeV gamma crisis', we examine a model whose one spatial 
direction is discretized at a high energy scale. Assuming the standard extra- 
galactic IR photon distribution, we evaluate the mean free-path of a energetic 
photon which acquires an effective mass in the model. For a wide range of 
the value of the lattice energy scale between a few TeV and around 10 10 GeV, 
the mean free-path of a TeV energy photon can be enlarged enough to solve 
the 'crisis'. Taking into account the effects of the universe expansion, we find 
bounds of the lattice constant. 
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1 Introduction 



Gamma-rays of TeV energies from distant BL Lacertae objects, or 'blazars', 
provide us a problem so-called 'TeV gamma crisis' lj. For example, HEGRA 
collaboration observed the gamma-ray from Markarian 501 (Mrk501) in 139 
Mpc distance, and the maximum measured energy of the gamma-ray reaches 
to about 20 TeV 0. An energetic gamma-ray can be attenuated by scatter- 
ings with the extra-galactic infrared photons in 77 — > e + e~ process. Since 
the evaluated mean free-path of a 20 TeV photon is only 13 Mpc, the gamma- 
rays should have been suppressed by a factor of exp(— 139/13) = 1/(4 x 10 4 ). 
It means that the original brightness of Mrk501 in 20 TeV is 4 x 10 4 times 
brighter than what HEGRA saw. Comparing the observed energy spectrum 
of the gamma-rays of Mrk501 with that of much closer TeV gamma-ray 
sources like the Crab nebula, such a strong suppression in high energy side 
in the energy spectrum seems to be unrealistic. The same problem also holds 
for Mrk421, 1ES1959+650 and H1426+428 0. The observed highest energies 
of the photons from these blazars are listed in ITablc II 

Some trials to solve the 'crisis' have been proposed by many authors. 
The most conservative way is to assume a less density of the extra-galactic 
infrared photons. One can also evaluate the appropriate IR photon density 
so that the traveling gamma-rays do not absorbed so much There is 
another possibility to assume the violation of the Lorentz invariance in the 
space-time 011]. Such kind of idea is motivated by a recent development of 
the theories of the large extra-dimension. In the deformed space-time, it is 
expected that the energy threshold of the 77 — > e + e~ process can be shifted 
to reduce the absorption of the TeV gamma-rays. 

In this paper, we examine the possibility to answer the 'TeV gamma crisis' 
by a model with reduced degree of space freedom at a high energy scale [H]. In 
a model whose one spatial direction is discretized at a high energy scale, we 
evaluate the extra-galactic absorption of an energetic photon which acquires 
an effective mass in the model. 

The paper is organized as follows. In the section |21 the adopted model is 
presented. The threshold shift of the scattering process and the mean free- 
path of the high energy photon in this model are evaluated in the section |3] 
The effects of the universe expansion on the photon propagation is discussed 
in the section HJ We give the conclusions in the last section. 
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2 The Model 



In the recent theories with large extra dimensions, we consider that the usual 
four dimensional space-time is located on a world volume of the brane (the 
extended object like membrane) and the world volumes stand parallel to 
each other with a separation distance (or lattice constant) of a. At high 
energy, such lattice structure becomes manifest itself and we are confined 
into the dynamics on a single brane, but after the universe cools down, the 
neighboring branes which are irrelevant so far, are possible to be dynamically 
connected and we arrive at the five dimensional theory at low energy. This 
was shown to happen in gauge theories in Following this idea, we are 
tempted to examine that the space-time dimensions become smaller and 
smaller according to the increase of energy scale in various theories as is 
discussed in jH]. 

Here, in order to solve the 'TeV gamma crisis', we look into the model 
in which the usual four dimensional QED works at low energy, but at high 
energy the lattice structure of one spatial dimension becomes manifest itself 
and we arrive at the three dimensional QED. More explicitly, we consider 
the theory in which the x 3 direction is discretized with the lattice constant 



a of the space structure. See Fig. 1 for an illustration of the model. 

If the discrete dimension is written as x 3 = n ■ a (n: integer), the Klein- 
Gordon equation for a boson field <fi becomes 

{(p-'lr) 2 -^ 2 }^^ 1 - 2 '") 

— \ U(t, x 1 , x 2 , n + 1) - 2(j)(t, x 1 , x 2 , n) + 0(t, x 1 , x 2 , n - 1)} ^ 
— m* 2 (f)(t, x 1 , x 2 , n) = 0, 

where t, x 1 and x 2 are the continuum time and space coordinates, m* is the 
boson mass. 

In the momentum space it becomes 

(E* -k\- k\ - k\ - m* 2 ) = 0, (2) 

where and (^1,^2,^3) are the energy and the momentum of the boson 
with, 

I 2 ■ a h /q\ 
k 3 = - sm . (6) 
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Now we can understand that /c 3 is replaced by k 3 for the discrete dimension. 
Such replacement is valid also for a photon. 

The fermion field ip describing electron satisfies the following Dirac equa- 
tion: 

f . d . i d -2^1,/ 1 2 \ 

+«7 3 ^ {ip(t, x 1 , x 2 , n + 1) - ip(t, x 1 , x 2 , n - 1)} ^ 
-Mip{t,x l ,x 2 ,n) = 0. 

In the momentum space it becomes 

( 7 °£ e + 7V1 + 7 2 P2 + 1% - M) V = 0, (5) 

where £7 e and {pi,P2,P3) are the energy and the momentum of the fermion 
with, 

p 3 = -sinap 3 . (6) 

As for the interaction of electron and photon, the gauge field (the con- 
nection field) A% in the x 3 direction is generated dynamically between the 
neighboring branes at x 3 = na and (n + l)a, giving 

ip(t, x , x 2 , 71)73 ?7 (n, n + l)if)(t, x 1 , x 2 , n + 1) + h.c, (7) 

with 

U(n, n + 1) = exp zeav4 3 (t, x 1 , x 2 , n) = 1 + ieaA 3 (t, x, x 2 , n) + • ■ ■ (8) 

In the expansion (jHJ), the first term contributes to the Dirac equation dis- 
cussed above, and the second 0(a) term contributes to the interaction. The 
lowest interaction reads 

iea7p(t, x , x 2 , n)7 3 A 3 (t, x 1 , x 2 , n)ip(t, x , x 2 , n + 1) + h.c. (9) 

The incoming fermion and the outgoing fermion are located on the different 
branes, so that the gauge interaction is modified from the usual QED by the 
following extra factor: 

- (e" iap3 + e aip a) , (10) 

where p' 3 and p 3 are, respectively, the third component of the incoming and 
outgoing fermion momenta. However, such modification of interactions from 
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the standard QED always associate the factor ((lattice constant a) x (mo- 
mentum or energy)) 2 , where the square guarantees the parity conservation. 
Therefore, the modification of the interaction is small at energy lower than 
a -1 . It is also expected that the QED in our model has additional vertices 
of photons with electrons. However, the effect of such deformation of the 
QED is expected to be less important for the 'crisis' problem than that of 
the shift of the threshold. The most important change of the QED for the 
'crisis' is the energy-momentum relation which provides the threshold shift 
of the 77 — > e + e~ process, as will be presented below in the section El In this 
paper, we concentrate on the effect of the change of the energy-momentum 
relation. 

In a way of thinking, particles are assumed to have always the usual 
energy-momentum relations such as, 



E k^ A; 2 



^3 



0. 



:n) 



But at high energy scale, the x 3 space is curled up to a circle with the circum- 
ference a. Then, in this deformed space, the equation of motion is written 
using &3 instead of k 3 in equation ([lip. Therefore, the energy-momentum 
relation can be viewed effectively as 



E^ k\ k\ k 3 



m 



*2 



;i2) 



having an effective mass 



m 



*2 



^3 



2 ak 3 
— — sin 

ak 3 2 



(13) 



A similar relationship also holds for a fermion with the intrinsic mass m e 



El 



2 2 
P2-P3 



7712 2 2-2 



m 



*2 



m, 

€ 

M. 



2 1 *2 

m e +m e , 



Pi 



ap 3 



sin ap 3 



(14) 
(15) 

(16) 



It is essential that the additional component of the effective mass m* only 
depends on a and p 3 , but does not depend on m e . Expanding the equations 
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(JT3J) and in powers of ak 3 and ap 3 , respectively, one gets, 

(17) 



m* 2 = kl 



*2 2 

m e = p 3 



— (ak 3 f + —(ak 3 f + - 
12 v J; 360 v ; 

^(ap 3 ) 2 + —(ap 3 ) 4 + ■ ■ 



(18) 



3 Threshold Effect to Mean Free-Path 



The scattering process of the blazar gamma-rays 77 — > e + e~ is characterized 
by an asymmetric initial photon energies; i.e., the energy of the traveling 
photon E 1 is in the order of TeV, while that of the extra-galactic IR photon 
e is in the order of eV. Thus, the effect of the deformation of the space 
structure to the extra-galactic IR photon is totally negligible. The energies 
of the produced electron E e \ and positron E e2 are also in the order of TeV 
in general. 

Especially, at the threshold, the electron and the positron have the same 
energy E e = E e \ = E e <i ~ E 1 /2. Assuming the blazar photon travels along 
to the x 3 direction, we can write down the energy-momentum conservation 
law just on the threshold as follows: 

£ 7 + e = 2E e , (19) 
E^-e = 2E e p e , (20) 



where /3 7 = y 1 — m* 2 /E^ and (3 e = \J 1 — M 2 / El are the velocities of the 
energetic photon and the electron in the laboratory frame, respectively. Since 
E 1 ~ k 3 and E e ~ p 3 , m* ~ 2m* at the threshold. In such a case, M e is 
dominated by m* for, 

a" 1 < J E 2 /(4 v / 3m e ), (21) 

which results in a significant increase of the threshold energy and a desirable 
suppression of the absorption process. For E 1 = 20 TeV, a -1 < 1 x 10 11 
GeV which means that the deformation of the space structure at a very high 
energy scale can affect the electron-positron pair production at the scale of 
MeV. 

The mean free-path of the energetic photon be evaluated by the 

equation (1) in the reference PQ: 



^-mL^i"**' - m * 2K (22) 
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where n(e) is the number density distribution of the extra-galactic IR photon 
of the energy between e and e + de which can also be found in the figure 1 in 
the reference [T], s — 2E 1 e{l + /3 7 ) + m* 2 is the CM. energy squared of the 
process, and, 

4M 2 -m* 2 ^ Ami m) 
2E 7 (l+^ 7 ) " 2E 7 (l+^ 7 )' 1 ] 

Smax = m* 2 + 2eE 1 (l + p), s min = 4M 2 , (24) 



where /? = 1 — m* 2 / s is the average initial photon velocity in the CM. frame. 
We used for simplicity the formula of the total cross section of the scattering 
process a as same as in the ordinary QED, which is described as, 



Ana 2 

a 



^4 



(25) 



where a the QED fine structure constant, /3 e = yl — 4M 2 /s the electron 
velocity in the CM. frame, respectively. We adopted an approximation such 
that the effective masses of the electron and the positron has the same value 
m* = m*/2 as just as at the threshold. We took into account the effective 
mass effect of the photon only in the flux factor in the denominator of the first 
fraction in the r.h.s. of the formula (|25|) . The cross section a is illustrated in 



Fig. 2 for the ordinary space-time with 0.511 MeV electron mass. The peak 
of the cross section can be found at 1.40 times above the threshold. 

|Fig. 3| shows the evaluated value of the mean free-path of the traveling 
photon x 77 as a function of the photon energy E~ and the lattice constant of 
the space structure a. The lowest curve of a -1 = oo case corresponds to the 
ordinary continuum space-time. Here the upper edges of the curves, which 
corresponds the lower bounds of a -1 for a fixed £7 7 , can be derived from a 
restriction of aE^/2 < ir, which proves the photon propagates to the correct 
direction. 

We found that the mean free-path is enlarged drastically for a wide range 
of the value of a -1 . One can observe from the figure, departure of the curves 
from the a" 1 = oo case happens around a -1 = E 2 j \A\f?>m e ) , as is already 
suggested in the equation (|2*Tj) . The mean free-path in the modified space 
structure is characterized by the followings: 

(i). A minimum of the mean free-path is found at E 1 around several hun- 
dred GeV to 1 TeV for a~ l < 10 9 GeV. 
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(ii) . Almost constant value of the mean free-path is found at around 10 

TeV and a -1 < 10 10 GeV, which is in contrast with the a -1 = oo case. 
It means that there is no energy cut-off of the blazar photon spectrum 
around 10 TeV. 

(iii) . For a large E y , the mean free-path is raised dramatically as is propor- 

tional to a ■ E^ 2 . 



The |Fig. 4 shows that a 20 TeV photon has a longer mean free-path than 



139 Mpc of Mrk501 distance for 10 3 - 5 GeV < a' 1 < 10 10 GeV. 



4 Photon Propagation in the Expanding Uni- 
verse 

More accurately, one can take into account the effect of the universe expan- 
sion during the photon propagation from Mrk501 to the earth. This effect is 
only important near the lower bound of a -1 which results in a larger value 
of m* and slower (3 T 

We assumed a de Sitter expanding universe which is spatially uniform at 
the temperature below a -1 . The scale factor R(t) of the universe has a time 
t dependence as, 

^ = H ^Q m /R + Q v R 2 , (26) 

where H Q the Hubble constant, Q m the present matter density fraction and 
Q v the present vacuum energy density fraction, respectively. We adopted H 
= 71 km/s/Mpc, Q m = 0.27 and Q v = 0.73, according to the recent results 
by WMAP jTUj. During the universe expansion, the energy of the traveling 
photon decreases due to the red-shift. We assumed that the evolution of the 
number density distribution of the extra-galactic IR photon is only due to 
the cosmological expansion. The lattice constant a is assumed to be constant 
during the time of interest of the photon propagation. 

We can evaluate the distances of blazar galaxies by the Hubble law, ne- 
glecting any proper motions in the co-moving frame of the universe. At the 
time when the visual light we see today on the earth was departed from the 
blazar at £*, the scale factor of the universe expansion was, 

R{U) = Rq/{1 + z), (27) 
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where z is the red-shift of the observed light in the visual wavelength, R is 
the present scale factor. The present distance of the blazar in the co-moving 
frame can be evaluated as, 

where we set the value of the present time to be zero. Since an energetic 
photon with an effective mass travels at a slower speed /3 7 than unity, such 
a photon should depart at an earlier time £* than t* of the ordinary light: 

* = jC*|iA (29) 

Note that j3 1 is a function of the scale factor R(t), then it depends on t. The 
degree of the transparency T of an energetic photon can be evaluated as, 



T = exp 

here x 77 is also a function of R(t). 
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(30) 



Fig. 4| displays the transparency T dependences on the lattice constant 



a. 



for a 20 TeV photon from Mrk501 (z = 0.0336) and a 5.37 TeV photon from 
H1426+428 (z = 0.129). 

For MrkSOl, we found 10 3 - 4 GeV = 2.5 TeV < a' 1 < 10 10 - 5 GeV. Here 
we set the upper bound of a -1 so that the transparency T should be at least 
around 1/10. To reach a 20 TeV photon from Mrk501 to the earth, aE 1 (t)/2 
should be greater than ir even at the time of departure when the photon 
energy E 1 has its maximum value. The lower bound can be derived from 
this condition. 

One should recall here that we have assumed that the energetic photon 
propagates to the third direction which is in the discretization. If the photon 
travels in the direction with the polar angle 6, the third component of the 
momentum has the fraction of cos 9 of the absolute value of the momentum, 
and thus the limits of a~ l described above should be multiplied by cos# in 
this case. 

The same argument can be repeated to the other blazars listed in lTablcTI 
Among them, a new bounds to a -1 can be derived from the 5.37 TeV photons 
from H1426+428, where we found 10 314 GeV = 1.4 TeV < a' 1 < 10 9 - 5 GeV. 
Note here that the angle between Mrk501 and 1H1426+428 is 27°, so the 
limits of two blazars cannot be combined into a single expression trivially. 



9 



5 Conclusions 



To solve the 'TeV gamma crisis', we have examined a model whose one spatial 
direction is discretized at a high energy scale. Assuming the standard extra- 
galactic IR photon distribution, we have evaluated the mean free-path of 
an energetic photon which acquires an effective mass in the model. For a 
wide range of the value of the lattice energy scale a -1 between a few TeV 
and around 10 10 GeV, the mean free-path of a TeV energy photon can be 
enlarged enough to solve the 'crisis'. Taking into account the effects of the 
universe expansion, we have found bounds of a -1 : 2.5 TeV < a -1 < 10 10,5 
GeV for a 20 TeV photon from Mrk501, and 1.4 TeV < a~ l < 10 9 5 GeV for 
a 5.37 TeV photon from H1426+428, under an assumption that the blazar 
gamma-ray travels to the discretized direction. 

Further test of the present model can be performed by watching PeV 
photons. In the ordinary space-time, a 1 PeV photon has the mean free- 
path of 9 kpc. However, it can travel over 300 Mpc in our model with a~ l 
< 10 10 GeV, as can be seen in Fig. 3 Observation of PeV photons from 



the near-by galaxies will support our model, which is a consequence of the 
item |(iii)1 described in the end at the sectional For a -1 ~ 10 8 -10 9 GeV, a 1 
TeV photon has shorter mean free-path than a few to several TeV photons 
in our model, as was also mentioned in the item |(i)| It results in a dent of 
the photon energy spectrum around 1 TeV for the objects which are distant 
as hundred Mpc. 

We also mention some other possible effect of the space discretization. In 
the case of smaller values of few TeV to several ten TeV, collider 

experiments may find some direct evidence that the energetic light particles 
fly slowly due to the effective mass. Modification to the QED vertex structure 
may affect to the magnetic and the electric dipole moments of the elementary 
charged particles, even in the interactions with low energy photons. These 
observations will bring us additional bounds of a -1 . 
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Table 



Table 1. The maximum energy of the gamma-rays from the blazars observed 
by HEGRA. We adopt the highest energy bin in the energy spectrum 
such that the differential energy distribution is measured well above its 
error. 



object 


red-shift 


maximum 


significance 


reference 


name 




photon energy 






Mrk501 


0.0336 


18.45 TeV 


3.3a 


121 


Mrk421 


0.031 


13 TeV (read from Figure) 


11 


1ES1959+650 


0.047 


10.94 TeV 


4.0a 





H1426+428 


0.129 


5.37 TeV 


4.0a 


U 



Figures 

Fig. 1. An illustration of our model with a discretized space direction £3. The 
size of the lattice spacing a is common for all gaps, and each membrane 
is numbered by an integer. 

Fig. 2. The collision energy yfs dependence of the total cross section of the 
process 77 — > e + e~ in the ordinary space-time. The peak of the cross 
section can be found at y/s = 1.43 MeV = 1.40 x(2m e ). 

Fig. 3. The mean free-path x 77 of the energetic photon as a function of the 
photon energy E 1 and the lattice constant of the space structure en- 
ergy scale a -1 . The lowest bold curve of a -1 = 00 corresponds to the 
ordinary continuum space-time. 

Fig. 4. The transparency T dependences on the energy scale of the lattice 
constant a -1 , for a 20 TeV photon from Mrk501 (solid curve) and a 
5.37 TeV photon from H1426+428 (dashed curve). The lower bounds 
of a -1 are plotted by a square (Mrk501) and a cross (H1426+428). 
We assumed, for each case, that the blazar is on the third direction of 
discretization. 
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Fig. 4 
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